Skipping orbits and enhanced resistivity 
in large-diameter InAs/GaSb antidot lattices 
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We investigated the magnetotransport properties of high mobility InAs/GaSb antidot lattices. In addition 
to the usual commensurability features at low magnetic fields we found a broad maximum of classical 
origin around 2.5 T. The latter can be ascribed to a class of rosetta type orbits encircling a single antidot. 
This is shown by both a simple transport calculation based on a 'classical' Kubo formula and an analysis 
of the Poincare surface of section at different magnetic field values. At low temperatures we observe weak 
l/B-periodic oscillations superimposed on the classical maximum. 
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Periodically modulated two-dimensional electron gases 
(2DEGs) offer the possibility to study electron motion 
in artificially tailored periodic potentials.El If the mod- 
ulation potential is strong enough to deplete parts of 
the 2DEG around the potential maxima the system is 
called an antidot (AD) array. Antidot lattices are model 
systems to study both the classical chaotic motion of 
electrons in the potentialJandscape as well as result- 
ing bandstructure effects.El Investigations so far have 
essentially focused on the low field regime where the 
low-temperature magnetoresistance reflects peaks when- 
ever the classical cyclotron diameter 2R C fits around a 
(geometry-dependent) specific number of antidots.H The 
last maximum appears in the resistivity when 2R C is 
equal to the period a of the array. Here we focus on novel 
structure appearing at slightly higher fields correspond- 
ing to B- values where 2R C < a holds. While antidot 
lattices are usually based on GaAs/AlGaAs heterojunc- 
tions, we used here the material system InAs/GaSb to 
fabricate antidot arrays. Due to the pinning of the Fermi 
energg within the conduction band at open surfaces of 
InAsj3 depletion regions around the antidots should be 
significantly smaller than in GaAs based systems. This 
should allow, due to the resulting steep potential posts, 
to fabricate very short period antidot arrays or, as it is 
the case here, the fabrication of large antidots with small 
constrictions between neighboring antidots. 

Our samples are fabricated from InAs/GaSb het- 
erostructures grown on undoped GaAs substrates by 
molecular beam epitaxy. The epitaxial layers on the sub- 
strate consist of a 1.1 /zm AlGaSb layer, a 0.5 /im thick 



GaSb buffer followed by a ten period AlSb/GaSb super- 
lattice, the GaSb barrier of 50 nm thickness, the 15 nm 
InAs quantum well and a 5 nm GaSb top layer (see upper 
inset in Fig. [l]). While the heterostructure is nominally 
undoped, the carrier concentration n s can be adjusted 
by the top layer thickness.Qtl Typical values of n s in our 
unpatterned samples were about 1.5 ■ 10 12 cm~ 2 . Hall- 
bars were definedJpy standard photolithography and wet 
chemical etchingE It turned out that the alkaline pho- 
toresist developer attacks the GaSb-layer and decreases 
the carrier mobility from about 400 000 cm 2 /Vs to about 
100 000 cm 2 /Vs. An annealing step at 270°C in a contin- 
uous flow of forming gas restored the mobility to about 
300 000cm 2 /Vs. This corresponds to a mean free path 
of 6 /xm which is much larger than the lattice period 
a = 380 nm, thus making commensurability features ob- 
servable. Since InAs does not form a Schottky barrier 
to metals, it was not necessary to use alloyed contacts. 
Instead of gold, which proved to react with GaSb in the 
final annealing step, we used evaporated indium contacts. 
The antidots were defined with electron beam lithogra- 
phy and etched with the solution described above. A 
scanning electron micrograph of the samples showed the 
antidots to be uniform in size and shape with an AD 
diameter of 250 nm. 

The four-point magnetotransport measurements were 
carried out at temperatures ranging from 1.5 K to 40 K 
using standard lock- in techniques. Utilising the negative 
persistent photo-conductivity of InAs quantum wells ,u 
the carrier density could be varied by illuminating the 
sample with a red LED. The magnetoresistance and Hall 
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resistance traces at a carrier density of 1.33 ■ 10 12 cm -2 
and a temperature of 1.5 K are shown in Fig. [l| At low 
fields two commensurability peaks in the magnetoresis- 
tance curve can be observed. The peak at about 1 T be- 
longs to trajectories going around one antidot, the other 
one at 0.25 T corresponds to a deformed cyclotron orbit 
around four antidots. In the simplest model, the en- 
hanced resistivity at these field values can be ascribed 
to electrons pinned on commensurate cyclotron orbits, 
which cannot carry current At fields above 3.5 T 
strong Shubnikov-de Haas (SdH) oscillations appear and 
from 6 T spin-splitting of the Landau levels is resolved. 
By comparison to an unpatterned reference sample, we 
can conclude that the quality of the 2DEG was not se- 
riously affected by the patterning process. In the region 
between 1.5 T and 3.5 T we find a new broad peak struc- 
ture whose origin will be discussed below. In earlier work, 
a corresponding shoulder was ascribed to etch-induced 
defectsjj For a sample with triangular antidots, a high- 
field shoulder at 2R C — a/2 was explained hjz. electrons 
being reflected on the straight antidot edgesO 

Figure || shows the dependence of the magnetoresis- 
tance traces on the carrier density (changed by illumina- 
tion) and temperature (inset). The novel peak structure 
shifts similar to the main commensurability maximum at 
2R C = a. This suggests, together with the weak temper- 
ature dependence of the peaks (see lower inset of Fig. ||), 
a classical origin of the phenomenon. Associating the 
magnetic field value at which the broad maxima start 
to appear (marked by arrows) with the matching condi- 
tion R c = b we find a value of b — 110 nm. This value 
b corresponds to the distance between neighboring anti- 
dots in our device. This finding is consistent with the 
assumption that the broad maximum is associated with 
rosetta shaped orbits skipping around one antidot, as we 
will show below. 

Superimposed on the classical peak we find 1/B- 
periodic quantum oscillations. While the periodicity of 
these oscillations is similar to the period of the SdH- 
oscillations at higher magnetic field, there is a distinct 
transition between the high-field oscillations and the ones 
superimposed upon the classical maximum. This can be 
seen in the upper inset of Fig. || where the filling fac- 
tor v is plotted versus l/B. At around 4T there is 
a pronounced kink in the otherwise linear slope which 
is observed at the transition between the classical peak 
regime and the quantum mechanical regime in which 
Shubnikov-de Haas oscillations dominate. Compared to 
the Shubnikov-de Haas oscillations the amplitude of os- 
cillations superimposed upon the classical peak structure 
is strongly suppressed. We speculate that the weak 1/B- 
periodic oscillations in the low field regime stem from 
quantization of the skipping orbits. The quantum oscilla- 
tions are suppressed at temperatures above a few Kelvin 
(see lower inset of Fig. |2|) while the classical features sur- 
vive up to temperatures of 40 K where they are blurred 
as the GaSb substrate becomes conducting. 

In order to clarify the origin of the classical maximum 
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FIG. 1. Magnetoresistance and Hall resistance of an 
InAs/GaSb antidot lattice. Lower inset: A SEM picture of 
the sample. Upper inset: Layer sequence of the heterostruc- 
ture as described in the text. Letters a to d denote B-field 
positions of interest, which are analysed numerically in Fig. ^. 

between 1.5 T and 3.5 T we carried out numerical calcu- 
lations following recent workJMa We evaluated the Kubo 
formulaE] 

e 2 m* f°° i 
irn Jo 

to calculate the conductivity tensor where the indices 
i and j stand for the x- and y-direction. The trajecto- 
ries and the resulting velocities Vi and Vj were calculated 
numerically assuming a hard wall potential. The brack- 
ets denote the average over phase space. Previously only 
chaotic trajectories-were considered in the evaluation of 
the Kubo formulaJlJ Since the features we are interested 
in appear at relatively high magnetic field values, most 
of the phase space is occupied by regular non-chaotic tra- 
jectories. Therefore we included both regular and chaotic 
trajectories in our calculations. By inverting the conduc- 
tivity tensor we obtain the experimentally determined 
diagonal resitivity p xx . The results of a simulation using 
the simple model with perfectly hard walls are shown in 
Fig. ||. The calculated traces are remarkably close to the 
experimental results: the additional—structure between 
1.5 and 3.5 T was nicely reproducedJlj 

To obtain information about the types of trajectories 
involved in the observed magnetoresistance anomalies 
we analysed the Poincare surfaces of section at different 
magnetic field values. This is done by injecting electrons 
at x — (between two antidots; see left hand side of the 
lower inset of Fig. |j) with velocity v y into the lattice and 
recording the v y and y values of the injected electrons 
intersecting the line x mod a = 0. To take into account 
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FIG. 2. Density dependence of the magneto-resistance. 
Electron densities (from top to bottom): 1.13 ■ 10 12 cm -2 , 
1.24 ■ 10 12 cm -2 , 1.33 ■ 10 12 cm -2 , and 1.42 ■ 10 12 cm -2 . Ar- 
rows indicate the beginning of the new structure, vertical lines 
mark the onset of strong SdH oscillations. Lower inset: Tem- 
perature dependence of the features at n s = 1.42 ■ 10 12 cm -2 . 
Temperatures are 1.5 K, 3K, 10 K, 20 K, and 40 K from bot- 
tom to top. Graphs are offset by 100 for clarity. Upper in- 
set: l/B-position of the minima of the quantum oscillations. 
Filling factors v correspond to the carrier density determined 
at high fields. A magnetic field axis is included for guidance. 



possible soft wall effects we assumed here an antidot po- 
tential of the form 



V (x, y) = Vq ( cos — • cos — 
V a a J 



with /3 — 16. The trajectories and velocities of the 
injected electrons were obtained by solving the classi- 
cal equations of motion numerically at several magnetic 
fields (indicated with letters a to d in Figs. |]andQ). 

At about IT (letter a), the commensurability condi- 
tion 2R C = a is satisfied, and the Poincare surface of 
section exhibits a stable island in Fig. ^a belonging to 
cyclotron orbits around one antidot. This leads to the 
fundamental commensurability peak in p xx . The mini- 
mum in the magnetoresistance at 1.9 T (letter b) is due 
to the fact that the large antidots leave no room for sta- 




magnetic field (T) 

FIG. 3. Magnetoresistance traces obtained with a simple 
simulation based on the Kubo formula. The commensurabil- 
ity condition 2R C = a is shown for n s = 1.3 ■ 10 12 cm -2 used 
in the calculation. Inset: A sketch of the hard-wall potential 
and a rosetta orbit encircling an antidot. Its magnetic field 
position is marked with an arrow in the large graph. 



ble trajectories. This minimum is not found in antidot 
lattices with smaller AD diameters. If the magnetic field 
gets higher, the Poincare surface of section c shows that a 
large fraction of phase space is occupied by periodic and 
quasi-periodic trajectories. These are rosetta shaped or- 
bits, which encircle one antidot and remain stationary 
even if an electric field is applied in the simulation, giv- 
ing rise to a maximum in the magnetoresistance at 2.5 T. 
The situation is quite similar to pinned orbits, which are 
responsible for the commensurability peak at 2R C = a. 
The high field behaviour (letter d, B = 4T) of the clas- 
sical structure is governed by the formation of unstable 
rosetta orbits. Instead of being composed more or less 
of semicircles, the small cyclotron diameter at high fields 
permits orbits which resemble circles revolving around 
one antidot. These orbits tend to drift perpendicular to 
an electric field, thus reducing p xx , even though the sta- 
ble orbits still occupy a large volume in phase space. 

When the magnetic field exceeds 3.5 T, the amplitude 
of the SdH-oscillations increases drastically. At lower 
fields the cyclotron diameter exceeds the width of the 
constriction between adjacent antidots. Therefore, elec- 
trons can move easily from one edge of the sample to 
the other by hopping from AD to AD. As soon as the 
cyclotron diameter becomes smaller, this backscattering 
process is no longer possible and the SdH-oscillations be- 
come much more pronounced. From the magnetic field 
position of this point we determine this critical cyclotron 
diameter to be 110 nm. If we assume a very steep anti- 
dot potential, this leads to an antidot diameter of 270 nm 
which is a little larger than the lithographic diameter of 
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FIG. 4. Poincare surface of sections at B-fields labelled 
from a to d in Fig. n[ y and v y are given in units of the 
lattice period a and v2uf, respectively. Trajectories shown 
below correspond to shaded areas in the Poincare surfaces 
of section. Note that for the highest magnetic field (d) the 
largest part of phase space is still occupied by stable rosetta 
orbits, only the orbits in the shaded region are unstable in 
an external electric field. Lower left corner: Position of the 
Poincare plane. 



250 rim. Since on InAs sprfaces the Fermi level is pinned 
in the conduction bandJ3 one would not expect a depleted 
region enlarging the antidots. However, the exact pin- 
ning position of the Fermi level depends on the detailed 
conditions of the etch-exposed surfaces and mayr-well be 
different from the bulk Fermi level in the 2DEGE-3 Con- 
sequently, the enhancement of the antidot diameter can 
be due to a softer potential. 

In summary, we fabricated and measured antidot ar- 
rays on InAs/GaSb with large AD diameters. We find a 
broad maximum in the magnetoresistance at intermedi- 
ate magnetic fields characterized by R c ~ b as lower and 
2R C ~ b as upper edge, which we show to be of classi- 
cal origin. The trajectories responsible for this maximum 
are identified to be rosetta orbits around single ADs. We 
also observe weak quantum oscillations on this maximum, 
which we speculate to originate from quantized rosetta 
orbits. Finally, we estimate the effective diameter of the 
ADs and find that it is comparable with the lithographic 



diameter. 
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